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Abstract: We analyze the supersymmetric contributions to the direct and mixing 
CP asymmetries and also to the branching ratios of the B — > Kir decays in a model 
independent way. We consider both gluino and chargino exchanges and emphasize 
that a large gluino contribution is essential for saturating the direct and mixing CP 
asymmetries. We also find that combined contributions from the penguin diagrams 
with chargino and gluino in the loop could lead to a possible solution for the branching 
ratios puzzle and account for the results of R c and R n within b — > 37 constraints. 
When all relevant constraints are satisfied, our result indicates that supersymmetry 
favors lower values of R c . Finally we study the correlations between the mixing 
CP asymmetry Sk°tt° and mixing CP asymmetries of the processes B — > <pK and 
B — > rj'K. We show that it is quite possible for gluino exchanges to accommodate 
the results of that observables 
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1. Introduction 

Recently the BaBar and Belle collaborations have measured the CP averaged branch- 
ing ratios and the CP violating asymmetries of B — > Kn decays [1-3]. These results, 
in addition to those from the B — > (pK and B — > rj'K, offer an interesting avenue to 
understand the CP violation and flavor mixing of the quark sector in the Standard 
Model (SM). 

In the SM, all CP violating observables should be explained by one complex phase 
5ckm i n the quark mixing matrix. The effect of this phase has been observed in kaon 
system. In order to account for the observed CP violation in this sector, 5ckm has to 
be of order one. With such a large value of 8ckm, the experimental results of the CP 
asymmetry of B —>■ J/tpKs are consistent with the SM. However, the experimental 
measurements of the CP asymmetries of B — > <pK, B — > rj'K and B — > Kir decays 
exhibit a possible discrepancy from the SM predictions. Furthermore, it is well known 
that the strength of the SM CP violation can not generate the observed size of the 
baryon asymmetry of the universe, and new source of CP violation beyond the Sckm 
is needed. 
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Decay channel 


BR xlO 6 


A CP 


Sf 


K°TT- 


24.1 ± 1.3 


-0.02 ±0.034 


— 


K-7T 


12.1 ±0.8 


0.04 ±0.04 


— 


K-7T+ 


18.2 ±0.8 


-0.113 ±0.019 


— 


K°7l° 


11.5 ± 1.0 


-0.09 ±0.14 


0.34 ±0.28 



Table 1: The current experimental results for the CP averaged branching ratios and CP 
asymmetries of B — * Kir decays. 



In supersymmetric extensions of the SM, there are additional sources of CP 
violating phases and flavor mixings. It is also established that the SUSY flavor de- 
pendent (off-diagonal) phases could be free from the stringent electric dipole moment 
(EDM) constraints [4] . These phases can easily provide an explanation for the above 
mentioned anomalies in the CP asymmetries of B — > <f>K and B — > rj'K [5-7]. We 
aim in this article to prove that in this class of SUSY models, it is also possible 
to accommodate the recent experimental results of B — > Ktt CP asymmetries and 
branching ratios. 

The latest experimental measurements for the four branching ratios and the 
four CP asymmetries of B — > Kit [1] are given in Table 1. As can be seen from 
this table, the measured value of the direct CP violation in B° — > K~ir + is A^JL + = 
— 0.113±0.019 which corresponds to a 4.2<r deviation from zero. While the measured 
value of A^+^o, which may also exhibit a large asymmetry, is quite small. As we will 
see in the next section, it is very difficult in the SM to get such different values for 
the CP asymmetries. 

Also from these results, one finds that the ratios R c , R n and R of B — > Kit 
decays are given by 



R r 



Rn — 

2 



BR(B + -> K + 7T°) + BR(B- -> 



_BR(B+ 
BR(B° 



K + n- 



BR(B- 
- BR(B° 



K°7T- 

K-n- 



R 



BR(B° -> K°tt°) + BR(B° 
BR(B° -> K + ir~) + BR(B° -i 
BR(B+ -> K°7t+) + BR{B- -> K°tt 



K-it+X 



1.00 ±0.08, (1.1) 



T. 



B 
T B o 



0.79 ±0.08, (1.2) 



0.82 ±0.06. (1.3) 



In the SM the R c and R n ratios are approximately equal, however, the ex- 
perimental results in Eqs. (|l.l|J1.2j ) indicate to 2.4a deviation from the SM predic- 
tion. On the other hand the quantity R is consistent with the SM value. Here 
t bI t b° — 1-089 ± 0.017. These inconsistencies between the A^ and the R c — R n 
measurements and the SM results are known as Ktt puzzles. 
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These puzzles have created a lot of interest and several research work have been 
done to explain the experimental data [8,9]. It is tempting to conclude that any 
new physics contributes to B — > Kn should include a large electroweak penguin in 
order to explain these discrepancies. In SUSY models, the Z penguin diagrams with 
chargino exchange in the loop contribute to the electroweak penguin significantly for 
a light right handed stop mass. Also the subdominant color suppressed electroweak 
penguin can be enhanced by the electromagnetic penguin with chargino in the loop. 
Therefore, the supersymmetric extension of the SM is an interesting candidate for 
explaining the Kn puzzles. 

It is worth mentioning that also new precision determinations of the branching 
ratios and CP asymmetries of B — > nn have been recently reported [2,3]. However, 
the SUSY contributions to B — > tttt, at the quark level, is due to the loop correction 
for the process b — > dqq, while the SUSY contribution to B — > Kn is due to the 
process b — > sqq. Therefore, these two contributions are in general independent and 
SUSY could have significant effect to B — > Kn and accommodates the new result, 
while its contribution to B — > nn remains small. Thus we will focus here only on 
SUSY contributions to B — > Kn. 

In this paper, we perform a detailed analysis of SUSY contributions to the CP 
asymmetries and the branching ratios of B — > Kn processes. We emphasize that 
chargino contribution has the potential to enhance the electroweak penguins and 
provides a natural solution to the above discrepancies. However, this contribution 
alone is not large enough to accommodate the experimental results and to solve 
the Kn puzzles. We argue that the gluino contribution plays an essential rule in 
explaining the recent measurements, specially the results of the CP asymmetries. 
Recall that other supersymmetric contributions like the neutralino and charged Higgs 
are generally small and can be neglected. The charged Higgs contributions are only 
relevant at a very large tan (3 and small charged Higgs mass. Therefore, we are going 
to concentrate on the chargino and gluino contributions only. 

The paper is organized as follows. In section 2 we study the CP asymmetries 
and the branching ratios of B — > Kn in the SM. We show that within the SM 
the Kn puzzles can not be resolved. In section 3 we analyze the supersymmetric 
contributions, namely the gluino and chargino contributions, to B — > Kn. We 
show that a small value of the right-handed stop mass and a large mixing between 
the second and the third generation in the up-squark mass matrix are required to 
enhance the chargino Z-penguin. Also a large value of tan (5 is necessary to increase 
the effect of the chargino electromagnetic penguin. 

Section 4 is devoted to the constraints on SUSY flavor structure from the branch- 
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ing ratio of b — > sj. New upper bounds on the relevant mass insertions are derived 
in case of dominant gluino or chargino contribution. A correlation between the mass 
insertions (<^r) 23 and (S^ L ) 32 is obtained when both gluino and chargino exchanges 
are assumed to contribute significantly. In section 5 the SUSY resolution for the 
R c — R n puzzle is considered. We show that it is very difficult to explain this puzzle 
with a single mass insertion contribution. We emphasize that with simultaneous con- 
tributions from gluino and chargino one may be able to explain these discrepancies. 

In section 6 we focus on the CP asymmetries in B — > Kn processes. We show 
that with a large gluino contribution it is quite natural to account for the recent 
experimental results of direct CP asymmetries. The SUSY contributions to the 
mixing CP asymmetry of B° — > K tt° is also discussed. Finally, section 7 contains 
our main conclusions. 



2. B — > Ktt in the Standard Model 

In this section we analyze the SM predictions for the CP asymmetries and the branch- 
ing ratios of B — > Ktt decays. The effective Hamiltonian of AB = 1 transition 
governing these processes can be expressed as 

if£ B=1 = -| \{CiQ\ + C 2 Ql + Y J ^ + C 7l Q 7l + C, g Q, g )+h.c. ) (2.1) 

v p=u,c i=3 

where \ p = V p bV* s and Cj are the Wilson coefficients and Qi are the relevant local 
operators which can be found in Ref.[10]. Within the SM, the b — > s transition can 
be generated through exchange of PU-boson. The Wilson coefficients which describes 
such a transition can be found in Ref. [10] 

The calculation of the decay amplitudes of B — > involves the evaluation of 
the hadronic matrix elements of the above operators in the effective Hamiltonian, 
which is the most uncertain part of this calculation. Adopting the QCD factorization 
[11], the matrix elements of the effective weak Hamiltonian can be written as 

(nK\H eff \B) = ^ E X ^K\ (T p + T™) \B), (2.2) 



where 



and 



p=u,c 



10 



(nK\T p \B) = J2^K)WQi\B)F, (2.3) 



i=i 

10 



(nK\T;™\B) = /WW. J>(ttX)- (2.4) 



1=1 
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The term 7p arises from the vertex corrections, penguin corrections and hard spec- 
tator scattering contributions which are involved in the parameters a^nK). The 
(ttK \Qi\B)p are the factorizable matrix elements, i.e. if any operator Q = j 1 j 2 , 
then (irK\Qi\B) F = (n^B) (K\j 2 \0) or {K^B) {it\] 2 \ty ■ The other term T p ann in- 
cludes the weak annihilation contributions which are absorbed in the parameters 
biijiK). Following the notation of Ref. [11] we write the decay amplitude of B — > Kn 
as: 



p=u,c 



Spu/32 + Oi\- ^ a l,EW + 03 + 03,EW 



(2.5) 



V2A B -^ n o K - = ^2 \Kk [fi P u(ui + fo) + a\ + Oi^EW + 03 + 03,ew] 



p=u,c 



+ X r A * 

p=u,c 
p=u,c 



A lf- 



$pu&2 + 2 a 3,EW 



8 pu oii + a p 4 + a p 4EW + 0[ 



l^3,EW 



(2.6) 
(2.7) 



p=u,c 



-ol\ + 



-a 



p 

4,EW 



-B p 



+ X ? A ^ 



p=u,c 



5pu^2 + 2 a 3,EW 



(2- 



Here the coefficients of the flavor operators a P (nK) and 0?(nK) are given in terms 
of the coefficients a P (nK) and h^(nK) respectively [11]. The parameter A n ^ (A^) 
is given by i^m^F^ 71 ^ f K ^. Note that the parameters 6j of the weak annihila- 
tion and hard scattering contributions contain infrared divergence which are usually 
parameterized as 



X 



A,H 



(1 



PA,He 



i<t>A, 



(2.9) 



where pa,h are free parameters to be of order one, 4>a,h £ [0, 2ir], and Ah = 0.5. As 
discussed in Ref. [5], the experimental measurements of the branching ratios impose 
upper bound on the parameter p A . If one does not assume fine tuning between the 
parameters p and <ft, the typical upper bound on p A is of order of pa ^ 2. 

Fixing the experimental and the SM parameters to their center values, one can 
determine the explicit dependence of the decay amplitudes of the B — > Kir on the 
corresponding Wilson coefficients. For instance, with 7 = 7r/3, and Pa,h and <Pa,h 
are of order one, the decay amplitude of B° — > K~n + is given by 

Abo^+k- x 10 8 ^ (1.05 - 0.02 i)d + (0.24 + 0.07 i)C 2 + (3.1 + 14.5 i)C 3 
+ (4.9 + 37.7 i)C 4 - (2.9 - 13.1 i)C 5 + (5.5 - 43.7 i)C 6 
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+ (1.7 + 10.4 i)C 7 + (5.8 + 36.5 i)C 8 + (2.8 + 12.7 i)C 9 (2.10) 
+ (0.6 + 35.5 i)C w - (0.0006 + 0.04 i)C$f - (0.04 + 2.5 i)C e 8 f / . 

Similar expression can be obtained for B° — > K°7r°: 

Abo^ko x 10 8 ~ (-0.14 + 0.3 i)C l + (0.4 + 0.2 i)C 2 - (2.2 + 10.5 i)C 3 

- (3.5 + 26.7 i)C 4 + (2.1 - 9.3 i)C 5 + (3.9 - 30.9 i)C 6 

- (1.7 + 37.02 i)C 7 - (1.9 + 1.7 i)C 8 + (1.3 + 46.6 i)C 9 (2.11) 
+ (2.2 + 28.8 i)C w - (0.0002 + 0.01 i)C^/ + (0.03 + 1.8 i)C$ f . 

The amplitude of B~ — > K~tc° can be written as 

Ab-^ok- x 10 8 ~ (0.9 + 0.06 i)Ci + (0.6 + 0.3 i)C 2 + (2.2 + 10.5 i)C 3 
+ (3.5 + 26.7 i)C 4 - (2.1 - 9.3 i)C 5 - (3.9 - 30.9 z)C 6 

- (2.7 + 32.4 z')C 7 - (5.4 - 17.8 i)C 8 + (1.9 + 51.6 i)C 9 (2.12) 
+ (2.4 + 42.8 i)C XQ - (0.0004 + 0.03 i)C^ f - (0.03 + 1.8 i)C^. 

Finally, the amplitude of B~ — > K°tt~ is given by 



Ag-^-jfo x 10 8 ~ (0.4 - 0.4 - (0.00004 - 0.02 z)C 2 + (3.1 + 14.9 i)C 3 
+ (4.9 + 37.7 i)C A - (2.9 - 13.1 i)C 5 + (5.5 - 43.7 i)C 6 

- (3.2 + 3.8 i)C r - (10.8 + 13.8 i)C 8 - (1.9 + 5.7 i)C 9 (2.13) 

- (0.3 + 19.7 z)Ci + (0.0003 + 0.02 i)C^ f - (0.04 + 2.5 i)Cf/ , 

where = CV 7 — |Cs — Cq and C%t = C% g + C5. The SM contributions to the 
Wilson coefficients of b — > s transition, which are the relevant ones for B —>■ Kir, are 
given by 



Cf w ~ 1.077, Cl M ~ -0.175, C| M ~ 0.012, C£ A * ~ -0.33, Cf 1 ~ 0.0095 

'10 



C 5M _ _ .039, Cf M ~ 0.0001, Cf M ~ 0.0004, Cf M ~ -0.01, Cf M ~ 0.0019, 



C£f ~ -0.315, Cf s M ~ -0.149. (2.14) 

From these values, it is clear that within the SM, the dominant contribution to 
the B —>■ Kn decay amplitudes comes from the QCD penguin operator Q 4 . However 
the QCD penguin preserves the isospin. Therefore, this contribution is the same for 
all the decay modes. Isospin violating contributions to the decay amplitudes arise 
from the current- current operators Q\ and QV; which are called 'tree' contribution 
and from the electroweak penguins which are suppressed by a power a/a s . As can 
be seen from the coefficients of CV_i in Eqs. ([2.10| - [2.13|) , the electroweak penguin 
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contributions to the amplitudes of B — > Kir could be in general sizable and non 
universal. However, due to the small values of the corresponding Wilson coefficients 
in the SM Q2.14|) , these contributions are quite suppressed. 



Note also that the Q\ contribution to A B o^ w + K - and A B -^o K - is one order of 
magnitude larger than its contribution to the other two decay amplitudes. Therefore, 
in the SM, the amplitudes A B o^+ K - an d Ab--*^k- can be approximated as func- 
tion of C\ and C4, while the amplitudes Ago_ f7r o.K-o and A B -^ n - K o are approximately 
given in terms of C4 only. It is worth noting that the difference between the coeffi- 
cients of C4 in the amplitudes A B o^ w + K - and A B -_> n o K - is just due to the factor of 
y/2 in Eq.(CT, which is the same difference between the corresponding coefficients 
in A B - _ >?!-- k° & n d — A B o^ n oxo. 

We are now in a position to determine the SM results for the CP asymmetries 
and the CP average branching ratios of B — * Ktc decays within the framework of the 
QCD factorization approximation. The direct CP violation may arise in the decay 
B — ► Kit from the interference between the tree and penguin diagrams. The direct 
CP asymmetry of B° — > K~tc + decay A^L^ is defined as 

a cp _ \A(B^K-n+)\ 2 -\A(B^K+n-)\ 2 



\A(B° -> K-n+)\ z + \A(B° K- 



7T 



and similar expressions for the asymmetries Ag% A^L^q and A^ . Also the 
branching ratio can be written in terms of the corresponding decay amplitude as 

BR(B -> Ktt) = ^^j\A(B - ^vr)| 2 -L, (2.16) 

OTT M B 1 tot 

where 

jp| = [(Ml - (m K + m^f) {Ml - (m K - m.) 2 )] 2 

2-Mf? 

The SM results are summarized in Tables 2 and 3. In Table 2, we present the 
predictions for the branching ratios of the four decay modes of B — > 7r.fr. We assume 
that 7 = 7r/3 and consider some representative values of pa,h and 4>a,h to check the 
corresponding uncertainty. Namely, pa,h = 0,1,3 and 4>a,h = 0(1) are considered. 
From these results, one can see that for pa,h £ [0, 1] the SM predicted values for 
the branching ratios of B — > Ktt are less sensitive to the hadronic parameters. 
Larger values of Pa,h enhance the branching ratios and eventually they exceed the 
experimental limits presented in table 1 for p > 2. It is also remarkable that the 
SM results for the BR{B~ -> K°n-) and BR(B° -> K°n°) are larger than the 
experiment measurements, while the results for BR(B~ — > K~tt ) and BR(B° — > 
K~n + ) are consistent with their experimental values. This discrepancy does not 
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Branching ratio 


Pa,h = 


PA,H = 1 & (f>A,H ~ 1 


p A ,H = 3 & (j)A,H ~ 1 


BR R07T - x 10 6 


31.06 


33.35 


43.92 


BR K -^ x 10 6 


17.31 


18.45 


23.36 


BR K - n+ x 10 6 


25.87 


27.98 


39.55 


BR ROn0 x 10 6 


11.41 


12.47 


18.66 


Rn 


1.13 


1.12 


1.059 


R c 


1.11 


1.106 


1.063 


Rn 


0.83 


0.838 


0.9 



Table 2: The SM predictions for the branching ratios of the four decay modes of B — > Kir 
with 7 = 7r/3 



seem to be resolved in the SM, even if we consider large hadronic uncertainties. 
The parameters R c and R n , defined in Eqs.( |l.l| , |1.2| ) as the ratio of the CP average 



branching ratios of B — > Kit exhibit this deviation from the SM prediction in a clear 
way. The results in Table 2 show that in the SM R c ~ R n > 1. However, the recent 
experimental measurements reported in Table 1, implies that R c ~ 1 and R n < 1. It 
is very difficult to have this situation within the SM. As emphasized above, in the 
SM the amplitudes of B — > Ktc can be approximately written as 

Abo^+k- - Oi + h i) d + (a a + & 2 i) <^4, (2.18) 
Aflo^jro ^ _ 7/|( a 2 + 6 2 i) C 4 , (2.19) 

A B -^o K - ~ -^(a! + 6i i) Cx + -J=(a 2 + 6 2 z) C 4 (2.20) 
Ab-^-^o ~ (a 2 + 6 2 i) C 4 . (2.21) 

Thus, the parameters i? c and _R n are given by 

R c = R n = 1(01 + 61 + + ^ ^ > 1 ( 2 . 22 ) 
|(a 2 + b 2 i) C4I 

which is consistent with the result given in Table 2, using the full set of the Wilson 
coefficients. 

Now we turn to the SM predictions for the CP asymmetries of B — > Kir. Let us 
start by considering the approximation that the decay amplitudes for B~ — > K°7T~ 
and B° — ► fC 7r° are dominated by the pure gluon penguin operator Q 4 while the 
amplitudes for B~ — > K~tt° and i?° — > K~n + are given by Q4 an d a l so by the tree 
contribution of the current- current operator Q\. In this case, the following results are 
expected: The direct CP asymmetries A^^- and A^o^o should be very tiny (equal 
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CP asymmetry 


Pa,h = 


p A ,H = 1 & <f>A,H ~ 


Pa,h = 3 & 4>a,h ~ 


aCP 


0.007 


0.0086 (0.005) 


0.0078 (0.001) 


aCP 
A K-k° 


0.029 


0.063 (-0.006) 


0.185 (-0.15) 


aCP 
/i K-iv+ 


0.0044 


0.057 (-0.049) 


0.194 (-0.19) 


aCP 
A RW 


-0.02 


-0.013 (-0.025) 


-0.019 (-0.002) 



Table 3: The SM predictions for the direct CP asymmetries of the four decay modes of 
B — > irK with 7 = 7r/3 



zero in the exact limit of this approximation). The direct CP asymmetries A^o^o 
and A^. n+ should be of the same order and larger than the other two asymmetries. 

The SM results of the CP asymmetries for the different decay modes, including 
the effect of all local operators Qi, are given in Table 3. As in the case of the 
branching ratios, we assume that 7 = 7r/3, and pa,h = 0, 1, 3. Respect to the strong 
phases </>a,h, we take it to be of order one as before. Due to the sensitivity of the CP 
asymmetry on their sign, we consider both cases of 4>a,h — Few comments 

on the results of the direct CP asymmetries given in Table 2 are in order: 

1. The CP asymmetries A^. and A^L + are sensitive to the sign <pA (note that 
<pH is irrelevant for these processes). On the contrary, the CP asymmetries 
Afto^- and A^o n o are insensitive to this sign. 

2. As expected, the results of the CP asymmetries A^ w - and A^o^o are very 
small even with large values of p A - 

3. The value of A C j[_^ and A C K P 7 . can be enhanced by considering large value 
of pa and one gets values for A^L + of order the experimental result given 
in Table 1. However, it is very important to note that in this case, the CP 
asymmetry A^-^ i s also enhanced in the same way and it becomes one order 
of magnitude larger than its experimental value. 

While a confirmation with more accurate experimental data is necessary, the 
above results of the branching ratio and the direct CP asymmetries of B — > irK 
show that within the SM the current experimental measurements listed in Table 1 
do not seems to be accommodated even if one considers large hadronic uncertainties. 
It is worth stressing that the QCD correction would not play an essential role in 
solving this Kn puzzles. Furthermore, since we are interested here in the ratio of the 
amplitudes, many of the theoretical uncertainties cancel. So it can not be the source 
of these discrepancies. 
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Another useful way of parameterizing the decay amplitudes can be obtained by 
factorizing the dominant penguin amplitude P, where P is defined as [12] 



Pe %Sp — a\ a 4:,EW + + 03,EW 



(2.23) 



In this case, one can write the above expressions for the decay amplitude as follows: 



KKkP [l + r A e i5 -e-^], 



V2A B -^ 0K - = \ C A^P [1 + (r A e i5A + r c e lS ) + r EW e iS ™] , 



A B°- 



K- 



-V2A 



KA^p 



1 + (r A e i5A + r T e i&T ) + 



EW t 



(2.24) 



l + (r A e idA + r T e l5T - r c e iSc ) e" i7 + r% w e i5 ™ - r EW e 



where 



r A e i5A = 

r T e i5T = 

r c e lSp = 

i&EW 



3 



(3 2 + «4 - -at t EW + #3 + $i,EW 



/P, 



3 3 



/P, 



a l + Rktt&2 + 7:(PKir a 3,EW + a 4,E\v) 



/P, 



rEwe 



r EW e iS ^ 



(RKn^EW + a 4,Ew) 



/P, 



EW 



3,EW) 



IP 



(2.25) 
(2.26) 
(2.27) 
(2.28) 
(2.29) 



Here we define A u /A c = e KM e * 7 , R Kn = A^/A^, and 5 A , S T , 5 C , Sew and 8% w as 
strong interaction phases. The SM contributions within the QCD facorization leads 
to the following results: 



(Pe iSp ) 



SM 



-O.lle 1 



0.05H 



(r A e t5A ) SM = 0.019e 



0.26i 



(r c e i5c ) SM = 0.186e 2 - 9 \ (r T e^) SM = 0.191e 2 - 9 \ 
(r EW e iSEW ) SM = 0.13e-°- 2 \ ( r ^ e ^-) SM = 0.012 e - 2 - 5i 



(2.30) 



As can be seen from this result, within the SM r A and r EW are much smaller than 
r c , r T and r EW , so that they can be easily neglected. In this case, the parameters 
R c and R n can be expressed by the following approximated expressions 



Rc - 



1 l-ra — 



1 + 2r c cos 8 C cos 7 + 2r EW cos 8 EW , 
1 + 2r T cos 5t cos 7 
1 + 2r T cos 8 T cos 7 — 2r c cos 8 C cos 7 — 2r EH / cos 5 E w ' 



(2.31) 
(2.32) 
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which confirms our previous conclusion that in the SM R n ~ R c > 1. Explicitly, 
using the results of Eq.( [2.30| ), one finds that 

R c = 1.08(1.45), R n = 1.13(1.6), R = 0.757(0.673) (2.33) 

for 7 = 7r/3(27r/3), which is quite close to the full result that we obtained in Table 

2, with pa ~ 1. 

Now, we would like to comment on the mixing CP asymmetry of B — > Kw. 
CP violation in the interference between mixing and decay can be observed as time 
dependent oscillation of the CP asymmetry. The amplitude of the oscillation in 
charmonium decay modes provides a theoretical clean determination of the parameter 
sin 2/3 of the unitary triangle. The SM predicts the 5-decay modes, dominated by a 
single penguin amplitude such that B — > <j)K, B — > rj'K and B — > K°tt° to have the 
same time dependent CP asymmetry equal to sin 2(3. Again this result contradicts 
the experimental measurement given in Table 1. Note that the latest experimental 
results on the mixing CP asymmetry of B — > <pKs process are given by [2,3] 

S^Ks = 0.50 ±0.25^1 (BaBar), 

= 0.06 ± 0.33 ± 0.09 (Belle) , (2.34) 

where the first errors are statistical and the second systematic. Thus, the average 
of this CP asymmetry is S^Ks — 0.34 ± 0.20. On the other hand, the most recent 
measured CP asymmetry in the B° — > rj'Ks decay is found by BaBar [2] and Belle 
[3] collaborations as 

SrfK s = 0.27 ± 0.14 ± 0.03 (BaBar) 

= 0.65 ±0.18 ±0.04 (Belle), (2.35) 

with an average S v >k s = 0.41 ± 0.11, which shows a 2.5<r discrepancy from the SM 
expectation. This difference among S^k, S v 'k, S k o w o and sin 2/3 is also considered 
as a hint for new physics beyond the SM, in particular for supersymmetry. 

3. B Ktt in SUSY models 

As mentioned in the previous section, due to the asymptotic freedom of QCD, the 
calculation of the hadronic decay amplitude of B — > Kit can be factorized by the 
product of long and short distance contributions. The short distance contributions, 
including the SUSY effects are contained in the Wilson coefficients Cj. 

The SUSY contributions to the b — > s transition could be dominated by the 
gluino or the chargino intermediated penguin diagrams [5]. It turns out that the 
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dominant effect in both contributions is given by chromomagnetic penguin {Qgg)- 
However in case of B —> Kit, it was observed that this process is more sensitive to 
the isospin violating interactions [8,9], namely the contributions from the electro- 
magnetic penguin (Q77) an d photon- and Z-penguins contributions to Q-j and Qg. 
Therefore, in our discussion we will focus only on these contributions, although in 
our numerical analysis we keep all the contributions of the gluino and chargino. 

For the gluino exchange, it turns out that the Z-penguin contributions to are 
quite small and can be neglected with respect to the photon-penguin contributions. 
At the first order in the mass insertion approximation, the gluino contributions to 
the Wilson coefficients CV 7) 8 g , Cj and Cg at SUSY scale Ms are given by 



Cl(M s ) = C 9 (M S 



2a*a 



9V2G F m1 3 



LlJ 23-^042 (^, X), 



C~UM S ) 



8a, 7r 



9y/2G 



777 " 

(5 d LL ) 23 M 3 (x) + (5L) 23 ^M 1 (x) 



(3.1) 
(3.2) 



(S d LL ) 23 [ -M 3 (x) + 3M 4 (x] 



m b V 3 



+ (S d LR ) 23 ^[^M 3 (x)+3M 2 (x] 



(3.3) 



where x = m^/m 2 q and the functions M\(x),M 2 (x) and Pijk(x,x) can be found in 
Ref. [13, 14]. The coefficients CV 7 ,8g and 67,9 are obtained from CV 7i89 and 67,9 re- 
spectively, by the chirality exchange L <-» R. As can be seen from Eqs.( [3.2| , |3.3|) , the 
term proportional to (Sf R ) 23 in the coefficients CV 7)8g has a large enhancement factor 
nig /m&. This enhancement factor is responsible for the dominant gluino effects in 
5-decays, although this mass insertion is strongly constrained from b — ► 57. Note 
also that, since the photon-penguin gives the same contributions to Cj and Cg, and 
we neglect the Z-penguin contributions, we have Cj = Cg. Finally, it is clear that 
the coefficients 67,9 is suppressed with respect to CV 7) 8 S by a factor a/in at least. 

It is worth mentioning that the mass insertion ( 5^)23 can be generated by the 
mass insertion (S'[ L ) 23 as follows [15] 



where 



Therefore, 



($LR)23 - (#Ll)23 (<51r)33, 



rrib^Ab — //tan/3) 



tan/3 ~ l(T 2 tan/3. 



m 



m 



&) 2 3^10- 2 tan/?(^ L ) 23 . 



Hence, for a moderate value of tan/5 and (<5^ L )23 ~ 0(0.1), one obtains (SlrJzz of 
order 10~ 2 , which can easily imply significant contributions for the S^k and also 
account for the different results between S^k and S v >k- Thus in our analysis we 
define 

feW = &)23 + &) 23 (^33- (3-4) 

It is important to stress that in case of (Sf jR )23 eS dominated by double mass 
insertions, we still call this scenario as LR contribution. This is due to the fact that 
the main SUSY contribution is still through the C$ g which is enhanced by the chirality 
flipped factor m g /mb. In the literatures [16], this contribution has been considered 
in analyzing the CP asymmetry of B — > <pK and it was called as LL contribution, 
as indication for the large mixing in the squark mass matrix and dominant effect of 
(5^)23. However, we prefer to work with the notation LR c r to be able to trace the 
effective operators that may lead to dominant contributions for different B decay 
channels. 

The dominant chargino contributions are found to be also due to the chromomag- 
netic penguin, magnetic penguin and Z-penguin diagrams. As emphasized in Ref. [5], 
these contributions depend on the up sector mass insertion (5^)32 and (5^)32 while 
the LR and RR contributions are suppressed by A 2 or A 3 , where A is the Cabibbo 
mixing. At the first order in the mass insertion approximation, the chargino contri- 
butions to the Wilson coefficients are given by [5] 

CnM s ) = ^-(4C x + D x ), (3.5) 

C!m = £(4(l-J^)C x + D x ), (3.6) 

Cl ( = M 7 , (3.7) 
C* = M g , (3.8) 

where the functions F = C x (Z-penguin), D x (photon-penguin), M 7 (magnetic- 
penguin), and M 9 (chromomagnetic penguin) are given by [5] 



(^32 + ^)31 R f L + (^l) 3 2 + A&) 3 i Y t R™. (3.9) 



The functions Rp L and R§ L , F depend on the SUSY parameters through the chargino 
masses (m Xi ), squark masses (m)and the entries of the chargino mass matrix. For the 
Z and magnetic (chromomagnetic) dipole penguins Rq L ' rl and R 1 ^^ respectively, 
we have [5] 

Rc L = E \Vn\ 2 P { c\xi) + E [UnVaU^P^ix^xj) 

i=l,2 i,j=l,2 
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+ \Va\ 2 \V j i\ 3 l--P§\x i ,x j )) , 

i=l,2 »j'=l,2 

+ V*V jl Pg\x h x j )), 

Rjtf-f.g = l^l| 2 x Wi PMi,g{ x i) ~ Yf> VilUi2 x Wi — ~Pm1 „{ X di 

j i 

= -^V^^^P^Ji,,^), (3.10) 



where Y& is the Yukawa coupling of bottom quark, x W i = m^r/m 2 , Xi = rn^Jrh 2 , 
Xi = m 2 /m 2 ., and Xu = m 2 c ./m~ . The loop functions Pj^^ can be found in 
Ref. [5]. Finally, U and V are the matrices that diagonalize chargino mass matrix. 

Notice that the terms in and R\i g which are enhanced by m Xi /mft in 

Eq.( |3.10| ) lead to the large effects of chargino contributions to CV 7 and C$ g , respec- 
tively. Also the dependence of these terms on Yukawa bottom enhance the LL 
contributions in C 77j8g at large tan/5. In the case of light stop-right, the function 
Rq of the Z-penguin contribution is largely enhanced. In order to understand the 
impact of the chargino contributions in B — > Kir process, it is very useful to present 
the explicit dependence of the Wilson coefficients C 7 $j 1)8g in terms of the relevant 
mass insertions. For gaugino mass M 2 = 200 GeV, squark masses rh = 500 GeV, 
light stop rhi = 150 GeV, /i = 400 GeV, and tan/3 = 10, we obtain 



C 7 X ~ 0.000002(^ L ) 32 - 0.000011(5^)31 - 0.000046(^ L ) 32 , (3.11) 
CI ~ 0.00000039(5£ L ) 32 + 0.000037(5^)31 + 0.000165(5^ L ) 32 , (3.12) 
C% ~ -0.011(5^)31 - 0.05(52 L )32 - 0.00043(5^ L )3i - 0.002(5^ L ) 32 , (3.13) 
C* ~ -0.0032(5^)31 - 0.0014(5^)32 - 0.0003(5^) 3 i - 0.0012(5^) 32 . (3.14) 



From these results, it is clear that the Wilson coefficient C*L seems to give the domi- 
nant contribution, specially through the LL mass insertion. However, one should be 
careful with this contribution since it is also the main contribution to the b — > 57, and 
stringent constraints on (6^)32 are usually obtained, specially with large tan/3. Fi- 
nally, as expected from Eq. Q3.10p , only LL contributions to Of and C$ g have strong 
depend on the value of tan /3. For instance with tan f3 = 40, these contributions are 
enhanced with a factor 4, while the result of C* 9 and LR part of C* 7 and C% g change 
from the previous ones by less than 2%. 
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4. On the constraints from BR{B — ► X s 7) 

In this section we revise the constraints on SUSY flavor structure which arise from 
the experimental measurements of the branching ratio of the B — ► X s7 [17]: 



2 x 1(T 4 < BR(b -> s 7 ) < 4.5 x 10~ 4 (at 95%C.L.). 



(4.1) 



In supersymmetric models, there are additional contributions to b — > 57 decay be- 
sides the SM diagrams with H^-gauge boson and an up quark in the loop. The SUSY 
particles running in the loop are: charged Higgs bosons (H^ or chargino with up 
quarks and gluino or neutralino with down squarks. The total amplitude for this 
decay is sum of all these contributions. As advocated in the introduction, the neu- 
tralino contributions are quite small and can be safely neglected. Also the charged 
Higgs contributions are only relevant at very large tan/3 and small charged Higgs 
mass. Therefore, we consider chargino and gluino contributions only to analyze the 
possible constraints on the mass insertions (^ B ) 32 and (6^)23, where A = L,R. 

Although the gluino contribution to b — > 57 is typically very small in models with 
minimal flavor structure, it is significantly enhanced in models with non minimal 
favor structure [18]. In this class of models, both chargino and gluino exchanges 
give large contribution to the amplitude of b — > s 7 decay, and hence, they have 
to be simultaneously considered in analyzing the constraints of the branching ratio 
BR(b -> 57). 

The relevant operators for this process are Q 2) Qi^-, and Q 8g . The contributions 
of the other operators in Eq.(2T) can be neglected. The branching ratio BR(b — > s 7 ), 
conventionally normalized to the semileptonic branching ratio BR exp (B —>■ X c eu) = 
(10.4 ± 0.4)% [19], is given by [20] 



BR (B -> X s7 ) = BR cxp (B -> X c ev) 
x (\D\ 2 + A) (1 + 5 np ) 



\v t * s v tb 


2 Qa 

uu era 






2 


ng(z)k(z) 



1 



8 a 

3~ 7T 



( m b) \ 



(4.2) 



with 



i=i 



/i 



D = cf) M + ?M ( + Y, CfO.) [r,W + 7,?' log ^ 

A = ( e - a .Mlo S i(7+21„ 1! S)/3. _ |C.(») M |2 + £lM £ CfV)Cf M/ijffl 



7T 



i<i=l 



where z = m c /m^ \i is the renormalization scale which is chosen of order mb, and 
p is photon energy resolution. The expressions for Cf \ C\ l \ and the anomalous 
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dimension matrix 7, together with the functions g{z), k(z), Ti{z) and fij(S), can 
be found in Ref. [20]. The term 5 np (of order a few percent) includes the non- 
perturbative l/m& [21] and l/m c [22] corrections. From the formula above we obtain 
the theoretical result for BR(B — * Xgj) in the SM which is given by 

BB NLO (B -> X s7 ) = (3.29 ± 0.33) x 10" 4 (4.3) 

where the main theoretical uncertainty comes from uncertainties in the SM input 
parameters, namely m t , a s (M z ), a em , rn c / m b, m b,Vij, and the small residual scale 
dependence. The central value in Eq.( [4.3[ ) corresponds to the following central values 
for the SM parameters m p t ole ~ mp(m z ) ~ 174 GeV, mf = 4.8 GeV, mP° le = 
1.3 GeV, /i = m b , a s {m z ) = 0.118, a~ 1 (m^) = 128, sin 2 W = 0.23 and a photon 
energy resolution corresponding to p = 0.9 is assumed. 

The SUSY contributions to the Wilson coefficients CV 7 and C%g at leading order 
are given in the previous section. In general, the SUSY effects in b — > 37 decay can 
be parameterized by introducing i? 7 8 and R 7j $ parameters defined at the electroweak 
scale as 

p _ (CV 7 ,8g - C^ig) 5 _ CV 7 ,8g ,s 
^7,8 - TTSM ' ^7,8 - -^SM~, l 4 - 4 J 

where C^ t sg include the total contribution while C^ 8g contains only the SM ones. 
Note that in Cj-y t sg, which are the corresponding Wilson coefficients for Q^^^g respec- 
tively, we have set to zero the SM contribution. Inserting these definitions into the 
BR(B — > Xgj) formula in Eq. ( |4.2j ) yields a general parametrization of the branching 
ratio [18, 23] 



BR(B -> X sl ) = BR SM (S -> X s7 ) (l + 0.681i?e (R 7 ) + 0.116 
+ 0.083i2e(i2s) + 0.025 \Re(R 7 R* 8 ) + Re(R 7 R* 8 



ji?7| 2 + |i?7| 



+ 0.0045 



I-Rsl +1-^8 



(4.5) 



From this parametrization, it is clear that C*7 7 would give the dominant new contri- 
bution (beyond the SM one) to the BR(B — > X s 7) . Using the allowed experimental 
range given in Eq. ( |4.1|) , one can impose stringent constraints on CV 7 , and hence on 
the corresponding mass insertions. It is also remarkable that R 7 and R 7 have differ- 
ent contributions to the BR(B —>■ X s j), therefore, the possible constraints on C*7 7 
and hence on the LL and LR mass insertions would be different from the constraints 
on C7 7 and hence on the RR and RL mass insertions, unlike what has been assumed 
in the literatures. Furthermore, since the leading contribution to the branching ratio 
is due to Re(Rj), the CP violating phase of CV 7 will play a crucial role in the possible 
constraints imposed by BR(B — > X s j). 
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X 


l&) 23 | 


\(4lU 


0.3 


(a) 0.0116 

(b) 0.0038 

(c) 0.0012 


0.0038 


1 


(a) 0.02 

(b) 0.006 

(c) 0.002 


0.006 


4 


(a) 0.006 

(b) 0.015 

(c) 0.0045 


0.016 



Table 4: Upper bounds of \{^'[r(rl)') 23 \ f rom b —> sj decay for rriq = 500 GeV and 



arg(# 



LR(RL) 



(23 



(a),7r/2 (b),7r (c) respectively. 



Note that the constraints obtained in Ref. [13], namely (5f jR )23 < 1-6 x 10~ 2 and 
(5f jL )23 is unconstrained are based on the assumption that the gluino amplitude is the 
dominant contribution to b — > sj, even dominant with respect to the SM amplitude. 
Although this a very acceptable assumption in order to derive a conservative con- 
straints on the relevant mass insertions, it is unrealistic and usually lead to unuseful 
constraint. The aim of this section is to provide a complete analysis of the b — > S7 
constraints by including the SM, chargino and gluino contributions. 

Let us start first with gluino contribution as the dominant SUSY effect to b —>■ S7 
decay. We assume that the average squark mass of order 500 GeV and we consider 
three representative values for x = (m,g/mq) 2 = 0.3, 1, and 4. We also assume that 
the SM value for BR(B — > X s / y) is given by 3.29 x 10~ 4 , which is the central value of 
the results in Eq.( |4.3j ). In these cases we find that both the mass insertions | 23 1 
and I (5^)23 1 are unconstrained by the branching ratio of b — > 57 for any values of 
their phases. The upper bounds on | (Sf jR )23\ and | (5^)23 1 from b — > 57 decay are 
give in Table As can be seen from these results, the limits on | ( 5^)23 1 are quite 
sensitive to the phase of this mass insertion, unlike the bounds on | {S RL )23 \- Also, 
as suggested by Eq. (|4.5| ), the bounds on LR coincides with the ones on RL only if 
arg(<5L R ( RL ))23 = 7r/2. Note that in this case Re(R-?) vanishes and the expression of 
the branching ratio is a symmetric under exchange R7 and R7. 

Now we consider the chargino contribution as the dominant SUSY effect to 
b — > 57 in order to analyze the bounds on the relevant mass insertions in the up 
squark sector. From the expression of in Eq. ( P-13[) , which provide the leading 
contribution to the branching ratio of b — >• 57, it is clear that one can derive strong 



-17- 



M 2 \m 


300 


500 


700 


900 




M 2 \m 


300 


500 


700 


900 




(a) 0.04 


0.065 


0.095 


0.14 






(a) 0.17 


0.28 


0.45 


0.65 


150 


(b) 0.14 

(c) 0.51 


0.24 
0.85 


0.37 


0.54 




150 


(b) 0.65 

(c) - 










(a) 0.053 


0.075 


0.1 


0.15 






(a) 0.24 


0.34 


0.48 


0.67 


250 


(b) 0.20 

(c) 0.70 


0.28 


0.4 


0.55 




250 


(b) 0.86 

(c) - 










(a) 0.07 


0.09 


0.12 


0.16 






(a) 0.32 


0.4 


0.52 


0.73 


350 


(b) 0.26 

(c) 0.92 


0.33 
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Table 5: Upper bounds of K^^^l (left) and |(^ L )3i| (right) from b — > 57 decay for 
tan/3 = 10 and {i = 400 GeV and arg(c>L L )32(3i) = ( a )> 7r /2 (b),vr (c) respectively. 



constraints on (5^)32 and (5^)31 and a much weaker constraints (essentially no 
constrain) on {5\ L )^2 and 31. The resulting bounds on 32 and (5^)31 as 
functions of the gaugino mass M 2 and the average squark mass fh are presented in 
Tables 5, for tan/5 = 10 and fi = 400 GeV. 

The results in Table 5 correspond to positive sign of [i. If one assumed nega- 
tive sign of [/,, the constraints on | 32 1 an d I(^ll)3i| with arg(5L R )32(3i) will be 
exchanged with the corresponding ones with arg(5L L )32(3i) + n. Thus, in Table 5, the 
results of case (a) will be replaced with the results of (c) and vice versa. For larger 
values of tan/3, the above constraints will be reduced by the factor (tan/3/10). Note 
also that, because of the 577(2) gauge invariance the soft scalar mass Mq is common 
for the up and down sectors. Therefore, one gets the following relations between the 
up and down mass insertions 

= \Y$km Wl) Vckm]^ ■ (4.6) 

Hence, 

VDa = (5l L ) 32 + 0(X 2 ) . (4.7) 

As a result, the constraints obtained from the chargino contribution to b — > S7 
transition on | (5^)32! can be conveyed to a constraint on |(<^ L ) 3 2| which equals to 
|(5^ L ) 2 3|, due to the hermiticity of (Mq) ll . This is the strongest constraint one may 
obtain on |(<5f L )23|, and therefore it should be taken into account in analyzing the 
LL part of the gluino contribution to the b — > s. 
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Figure 1: Contour plot for BR(b — > 57) x 10 4 as function of (Sf jR )23 and (5^)32- 

Finally we consider the scenario in which both gluino and chargino exchanges 
are assumed to contribute to b — > S7 simultaneously with relevant mass insertions, 
namely (5f R )23 and (5^)32- It is known that these two contributions could give 
rise to a substantial destructive or constructive interference with the SM amplitude, 
depending on the relative sign of these amplitudes. Recall that in minimal super- 
symmetric standard model with the universality assumptions, the gluino amplitude 
is negligible, since (Sf R )23 < O(10~ 6 ), and the chargino contribution at large tan/5 is 
the only relevant SUSY contribution. In this class of model, depending on the sign 
of n the chargino contribution gives destructive interference with the SM result. 

In generic SUSY model, the situation is different and the experimental results of 
the branching ratio of b — * 37 can be easily accommodated by any one of these con- 
tributions. Also since the gluino and the chargino contributions are given in terms 
of the parameters of the up and down squark sectors, they are, in principle, inde- 
pendent and could have destructive interference between themselves or with the SM 
contribution. We stress that we are not interested in any fine tuning region of the 
parameter space that may lead to a large cancelation. We are rather considering the 
general scenario with large down and up mass insertions favored by the CP asymme- 
tries of different B processes. In this case, both gluino and chargino contributions 
to b — > S7 are large and cancelation of order 20 — 50% can take place. 
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Now, it is clear that the previous constrained obtained on (5 d R ) 23 and (S1 L ) 32 in 
Tables 4 and 5 will be relaxed. We plot the corresponding results for the correlations 
between (S d R ) 23 and (<^ L ) 32 in Fig. 1. Here we consider the relation (<^ L ) 23 = (^ L ) 32 
into account and also set (S RL ) 32 to zero. The phases of (S d R ) 23 and (S^ L ) 32 are 
assumed to be of order 7r/2 as favored by the CP asymmetry of B — > <pK. From this 
plot, we can see that constraints on these mass insertions, particularly (S^ L ) 32 are 
relaxed. 

5. SUSY solution to the R c — R n Puzzle 

Now we analyze the supersymmetric contributions to the B — > Kn branching ra- 
tio. We will show that the simultaneous contributions from penguin diagrams with 
chargino and gluino in the loop could lead to a possible solution to the R c — R n 
puzzle. As mentioned in section 3, these penguin contributions have three possible 
sources of large SUSY contribution to B — > Kir processes: 

1. Gluino mass enhanced Oy 7 and 8g which depend on (5^)23 and (S RL ) 23 . 

2. Chargino mass enhanced and 8g which depend on tan / 9(5^ L ) 23 . 

3. Right handed stop mass enhanced Z penguin which is given in terms of (S RL ) 32 . 

For the same inputs of SUSY parameters that we used above: m g = 500 GeV, 
m~ q = 500 GeV, m~ tR = 150 GeV, M 2 = 200 GeV, li = 400 GeV, and tan/3 = 10, one 
finds the following SUSY contributions to the amplitudes of B — > Kir 

A B o^ ORO x 10 7 ~ -9.82 i [(5 d LR ) 23 + (5 d RL ) 23 \ + 0.036 i (5 U LL ) 32 - 0.02 i (6 U RL ) 32 , 
Abo^c+k- x 10 7 ^ 14.04 1 [(5 d LR ) 23 + (4 L ) 23 ] + 0.06 % (5 U LL ) 32 - 0.001 1 (5 U RL ) 32 , 
Ab-^k- x 10 7 ~ 9.9 i [(S d LR ) 23 + (S d RL ) 23 \ - 0.04 i (S U LL ) 32 + 0.024 i {5 U RL ) 32 , 
A B -^- K o x 10 7 ~ 13.89 % [(5 d LR ) 23 + (5 d RL ) 23 ] + 0.05 1 (5 U LL ) 32 - 0.006 % (5 U RL ) 32 . 

It is remarkable that for the amplitudes Ago^^o^-o and Ab-^^k--, which suffer 
from a large discrepancy between their SM values and their experimental measure- 
ments, the SUSY contributions have the following features: (i) the effect of (5 RL ) 32 is 
not negligible as in the other amplitudes, (ii) there can be a distractive interference 
between the (5^)23 and (S1 L ) 32 contributions. As we will see below, these two points 
are important in saturating the experimental results by supersymmetry. Also note 
that the effect of gluino contribution through 0-j 1 is very small and the contribution 
of (Sf R ) 23 is mainly due to 8g . However the chargino effect of 7l can be enhanced 
by tan j3. 
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Figure 2: R c -R n correlation in SUSY models with |(<^) 32 | ^ 1, 1(^)23 1 6 [0.001, 0.01] 
and | (5^) 32 1 £ [0.1, 1]; see the text for the other parameters. The small and large ellipses 
correspond to la and 2a experimental results, respectively. 



We present our numerical results for the correlation between the total contri- 
butions (SM+SUSY) to the R n and R c in Fig. We have scanned over the rele- 
vant mass insertions: (o LL )32, (^lr)^ an d (^rl)32, since we have assumed (o LL )32 — 
(5 d LL ) 23 and (6 d LR ) 23 ~ (4 L ) 23 . We considered |(^ L ) 32 | G [0.1,1], |fe) 23 | G [0.001,0.01], 
arg[(<5£ L ) 32 ] e [-7T, vr],arg[(^ i? ) 23 ] ~ tt/3 (which is preferred by S^), and (<^ L ) 32 ] = 

1 (in order to maximize the difference between R n and R c ). As can be seen from 
the results in Fig. |2|, the experimental results of R n and R c at 2a can be naturally 
accommodated by the SUSY contributions. However, the results at la can be only 
obtained by a smaller region of parameter space. In fact, the values of R c is predicted 
to be less than one for the most of the parameter space. Therefore, it will be nice 
accordance with SUSY results if the experimental result of R c goes down. 

In order to understand the results in Fig. [| and the impact of the SUSY on the 
correlation between R n and R c , we extend the parametrization introduced in section 

2 for the relevant amplitudes by including the SUSY contribution [8]. In this case, 
Eqs. fl2.24|) can be written as 
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The parameters 5 a, 5c, $t, 5ew, 5 ew an d 8p, @ew, @ew are the CP conserving (strong) 
and the CP violating phase, respectively. Note that the parameters P,r EW ,r EW are 
now defined as 

Pe l9p e %5p = a\ — -oi% EW + (3% + (3^ EW , 



r EW e i(,EW e i5EW 



3 

2 {RKirOL3, E W + a 4,Ew) 



r EW e ie ^ w e iS ^ w 



o { a A,EW @l,Ew) 



/P, 

/P. (5.5) 



3 
2 1 

First let us include some assumptions to simplify our formulae. As is mentioned 
before, c^,a^ EW ,a^ EW ,P^,P^ EW ,Pl EW receive SUSY contributions through the 
Wilson coefficients. The upper index p takes both u and c, however the contribution 
with u index is always suppressed by the factor e KM — 0.018 so that its SUSY 
contributions can be safely neglected comparing to the one with the index c. As a 
result, (r^e* 5 - 4 ), (r c e t5c ) and (r T e l&T ) receive a correction of a factor 1/|1 + \ ■ 

Secondly we assume that the strong phase for SM and SUSY are the same. 
We found that this is a reasonable assumption in QCD approximation in which the 
main source of the strong phase comes from hard spectator and weak annihilation 
diagrams. This leads us to the following parametrization: 

Pe i5p e Wp = P SM e i5p {l + ke i6 'r) (5.6) 

r EW e i&EW e idEW = (r EW f M e i&EW (1 + le ie '™) (5.7) 

r^ w e i5 ^ e ^ C Ew = ( r c EW ) SM e iS ™ (1 + m/™). (5.8) 

where 

, ie' p _ ( a t ~ \ a< i,EW + ffis + ^,ew)svsy , , 

— / c_lc I qc I qc \ ' K^") 

\Ot A 2 a 4,EW ' P3 ' P3,Ew)sm 

l e i$' q _ (RkttCX^eW + q 4,.Ew)susy ^ ^ 

(RKTT a 3 >E W + a 4,£w)sM 

ifl' ( a l,EW ~ P3,E\v)svsy . . 

me = ~t~^ — ^ — (5-11) 

\ a A,EW ~ P3,EW)sm 

The index SM (SUSY) mean to keep only SM (SUSY) Wilson coefficients in a$ {>EW) 
and Pi( EW y Using these parameters, we also have 

rA JSA _ ( r Ae l5A ) S M iSc _ (r c e lSc ) SM iSr _ (r T e l5T ) SM 

TaC ~ h , ; i0' I' c _ h , ; id' I ' T ~~ h , ; iff \ \°- LZ ) 

1 1 + ke w p I 1 1 + fce 1 ^ I 1 1 + ke w p 

Now let us investigate the R c —R n puzzle. We shall follow the standard procedure; 
to simplify and expand the formulae. Considering the numbers obtained above, we 
shall simplify our formulae by assuming 
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1. the strong phases are negligible, i.e., 5p, 5 a, So, Sew, ^ew are an zero - 

2. the annihilation tree contribution is negligible, i.e. ta — 

3. the color suppressed tree contribution is negligible, i.e. rce t5c ~ rpt 1 ^ . 

Using these assumptions, we expand R c , R n and R c — R n . We expand in terms of rp 
and tew and r EW up to the second order. As a result, we obtain 

R c ~ l + r^-2r T cos('y + 9p) + 2rEwcos(6p-6Ew)-^rTrEwcos('y + 6Ew) (5.13) 
R c -R n ~ 2r T r EW cos( 7 + 26 P - 6 EW ) - 2r T r c EW cos( 7 + 26 EW ~ EW ) ( 5 - 14 ) 

Now, let us find the configuration which lead to R c — R n > 0.2. Looking at Eq. 
( |5.14p , we can find that in general, the larger the values of r T , r E w an d r EW are, the 



larger the splitting between R c and R n we would acquire. The phase combinations 
6p — 6 E w an d Op + 7 also play an important role. The possible solution of R c — R n 
puzzle by enhancing rpw, which we parameterize as I, has been intensively studied 
in the literature [9]. As we will see in the following, rp can also be enhanced due 
to the factor ke %e ' p which contributes destructively against the SM and diminish P. 
However, since P is the dominant contribution to the B — > Kir process, the branching 
ratio is very sensitive to ke %e ' p . Therefore, we are allowed to vary ke l6 ' p only in a range 
of the theoretical uncertainty of QCD factorization, which gives about right sizes of 
the B — > Kir branching ratios. As showed in Ref. [8], we would be able to reduce 
P at most by 30 %, which can be easily compensated by the error in the transition 
form factor p B ^-. K . 

Considering the tiny effect from the second term in Eq. (|5.14j), in order to achieve 



R c — R n > 0.2, we need rpTEw larger than about 0.1 or equivalently, Tew larger than 
about 0.5 with r^ M . In Ref. [8], it was emphasized that with k = 0, one needs / > 2 to 
reproduce the experimental values while an inclusion of a small amount of k lowers 
this bound significantly. For the SUSY parameters that we have considered above, 
the following results for our SUSY parameters k, I, and m are obtained 

ke i6p = -0.0019 tan (3(5 U LL ) 32 - 35.0(6 d LR ) 23 + 0.061(^) 32 (5.15) 
le ie " = 0.0528 tan P{5 U LL ) 32 - 2.78(6 d LR ) 23 + 1.11(^) 3 2 (5.16) 
me ie "o = 0.134 tan/3(^ L ) 32 + 26.4(^) 23 + 1.62(^) 32 (5.17) 

Note that we do not consider (6 d 3 )rl here but it is the same as (Sf R ) 23 with an 
opposite sign (see also [7]). Let us first discuss the contributions from a single mass 
insertion (S^ L ) 32 , {5j jR ) 23 or (S RR ) 32 to {k,l,m}; keeping only one mass insertion 
and switching off the other two. In this case, one finds that the maximum value of 
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{k, I, m} with 1(5^)32! = 1 is {k, I, m} = {0.061, 1.11, 1.62}. Thus, in this case where 
k is almost negligible, we would need I ~ 2 to explain the experimental data. We 
have a chance to enlarge the coefficients for (S^ R )s2 by, for instance, increasing the 
averaged squark mass rhq. However, even if we choose rhq = 5 TeV, we find that I is 
increased only by 20 to 30 %. The maximum contributions from {Sf jR )23 and {6^1)32 
are found to be {k, I, m} = {0.18, 0.014, 0.13} and {0.0019, 0.053, 0.13}, which are far 
too small to explain the experimental data. The coefficients for {Sf R )23 depend on 
the overall factor 1/rhq and on also the variable of the loop function x = m g /rhq and 
we found that m- g = rhq = 250 GeV can lead to 100 % increase. However, the value 
of I is still too small to deviate R c — R n significantly. As a whole, it is extremely 
difficult to have R c — R n > 0.2 from a single mass insertion contribution. 

Let us try to combine two main contributions, (5^)23 and (5^)32 terms. Using 
the previous input parameters and including the b — ■> 57 constraint |(<^r)23|, the 
maximum value is found to be {k,l,m} = {0.24, 1.12, 1.48}. In this case, it is easy 
to check that the experimental data are not reproduced very well [8]. As discussed 
above, for a large value of the averaged squark masses, I increases while k decreases. 
On the contrary, k also depends on the ratio of gluino and squark masses. Hence we 
need to optimize these masses so as to increase k and / simultaneously. For instance, 
with m~g = 250 GeV and rh- q = \ TeV, we obtain {k, /, m} = {0.30, 1.36, 1.90} which 
leads to a result within the experimental bounds of R c and R n . Finally, we consider 
the case with the three non-zero mass insertions. The main feature of this scenario 
is that we expect a relaxation of the constraints on |tan/3 x (<52l)32| and | (5^)23! 
from the cancelation between (5^)23 and (5^)32 contributions to b — > 57. Under 
this circumstance, we observe much larger R c — R n for various combination of the 
phases in this scenario. 

6. SUSY contributions to the CP asymmetry of B — > Ktt 

We start this section by summarizing our convention for CP asymmetry in B — > Kn 
processes. The time dependent CP asymmetry for B — > Kit can be described by 



where A Kn and Sk-k represent the direct and the mixing CP asymmetry respectively 
and they are given by 



A Kn (t) = A K7T cos(AM Bd t) + S K7r sm(AM Bd )t, 



(6.1) 



A 



Ktt — 



\p(Kn)\ 2 -l 
\p(K7i)\ 2 + V 



2Im{p(Kir)) 




(6.2) 
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where p(Kn) = e~ l ^ B ■ The phase (f>B is the phase of M12, the B° — B° mixing 
amplitude. The A(Ktt) and A(Kir) are the decay amplitudes for B° and B° to Ktx, 
respectively. 

The SM predicts that the direct and mixing asymmetry of B — > X7T decay are 
given by 

S K7r = sin2/3, C^ = 0. (6.3) 

The recent measurements of the CP asymmetries in B — > i^7r, reported in Table 
1, show significant discrepancies with the SM predictions. As mentioned above, 
SUSY can affect the results of the CP asymmetries in B decay, due to the new 
source of CP violating phases in the corresponding amplitude. Therefore, deviation 
on CP asymmetries from the SM expectations can be sizeable, depending on the 
relative magnitude of the SM and SUSY amplitudes. In this respect, SUSY models 
with non-minimal flavor structure and new CP violating phases in the squark mass 
matrices, can generate large deviations in the B — > Kit asymmetry. In this section 
we present and discuss our results for SUSY contributions to the direct and mixing 
CP asymmetries in B — > Kit. 

6.1 SUSY contributions to the direct CP asymmetry in B — > Kix 

Using the general parametrization of the decay amplitudes of B — > Ktt given in 
one can write the direct CP asymmetries A^? as follows: 



A C K P .^+ ~ 2r T sin S T sm(6 P + 7) + 1r c EW sin b G EW sin(# P - 9 C EW ) - r| sin 2S T sin 2{9 P + 7) 
+ 2r T r% w sin(S EW - S T ) sxn(9 EW + 7) - Ar T r% w sin 5 EW sm(9 P - 9 EW ) cos 5 T 
cos(#p + 7) - ^r T r% w sin S T sin(9 P + 7) cos 5 C EW cos(9 P - 9 EW ), (6.4) 
A° K ^- ~ 2r A sin 5 A sin(9 P + 7), (6.5) 
A KW — 2r EW sin S ew sin(6'p - 9% w ) - 2r EW sin 5 EW sm(9 P - 9 EW ), (6.6) 

— sin 5 T sin(9 P + 7) — 2r E w sin 5 EW sin(9 P — 0#vr) — r T sin 2£ T sin 2(6> P + 7) 

- 2r T r EW sm(5 EW - 5 T ) sin(9 EW + 7) - Ar T r EW sin 5 W sin(# P - 9 EW ) cos <5 T 
cos(0p + 7) — 4r T r B vy sin<5r sin(6 l p + 7) cos cos(9 p — 9 EW ). (6.7) 

From these expressions, it is clear that if we ignore the strong phases, then the direct 
CP asymmetries would vanish. However, Belle and BaBar collaborations observed 
non-zero values for the i^, thus we should consider no n- vanishing strong phases 
in this analysis. It is also remarkable that the leading contributions to the direct 
CP asymmetries are given by the linear terms of = rj-, ta, r E Wi r EWi un hke the 
difference R c — R n which receives corrections of order r^j. As in the previous section, 
we have assumed that the color suppressed contributions are negligible i.e., rce l5c = 
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rpe lST and we have neglected terms of order rf except for tt which is typically larger 
than r EW , r% w , and r A . 

The rescattering effects parameterized by ta are quite small [r s A M ~ 0(0.01)) 
therefore the CP asymmetry in the decays B — > K ^ is expected to be very small 
as can be easily seen from Eq. Qb.5| ). This result is consistent with the experimental 
measurements reported in Table 1. The sign of this asymmetry will depend on the 
relative sign of sin 5 a and sin(0p + 7). Note that the value of the angle 7 is fixed 
by the CP asymmetry in B — > tttt to be of order tt/3. The angle Op can also be 
determined from the CP asymmetry S^^k- 

In the SM, the parameters Ta,t ew are much smaller than r T ,r EW and Op = 0, 
therefore the following relation among the direct CP asymmetries Aj^ is obtained 

a cp > a cp > a cp . a cp 

This relation is in agreement with the numerical results listed in Table 3 for the direct 
CP asymmetries in the SM with pa,h,4>a,h — 1. To change this relation among the 
CP asymmetries and to get consistent correlations with experimental measurements, 
one should enhance the electroweak penguin contributions to B° — > K~n + decay 
amplitude, parameterized by r EW . Furthermore, a non- vanishing value of Op, which 
is also required to account for the recent measurements of S^k s an d S v >k s , is favored 
in order to obtain A^L^ > - It is worth mentioning that in the SM and due 
to the fact that Op = the second term in Eq. ( |6.4|) and Eq. (|6.7|) give destructive 
and constructive interferences respectively with first terms. Thus one finds A^^L^ is 
larger than A^. + . In SUSY models, the gluino contribution leads to a large value 
of Op and depending on the sign of this angle the parameter rp could be enhanced 
or reduced, see Eq. (|5.12 ). As will be seen below, in this case we can explain the CP 



asymmetry results with moderate values of the electroweak penguin parameter r EW . 
Note that in other models studied in the literatures, the value of this parameter is 
required to be larger than one in order to account for the CP asymmetry results. 

Now let us discuss the SUSY contribution to the CP asymmetries A^.. As can 
be seen from Table 1 that the experimental measurements of suffer from a 
large uncertainty. It turns out that it is very easy to have the SUSY results for this 
asymmetry within the range of 2a measurements. Thus, this decay mode is not useful 
in constraining the SUSY parameter space and can be ignored in our discussion for 
the correlation among the CP asymmetries of B — > Kit in generic SUSY models. 

We will consider, as in the previous section, three scenarios with a single mass 
insertion, two mass insertions, and three mass insertions. In the first case, if we 
consider the contribution due to the mass insertion (^2^)32 the maximum values of 
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{k, l,m} are given by {0.061,1.11,1.62}. While from (6^)23 and (S LL )32 one finds 
that the maximum values of {k,l,m} are {0.18,0.014,0.13} and {0.0019,0.053,0.13} 
respectively. Note that k is almost negligible in the case of dominant chargino con- 
tribution which depends on (6^)32 and ( 5^)32 an d can be significantly enhanced by 
the gluino contribution that depends on {S LR )23 as emphasized in Ref.[5]. Also from 
Eqs.(|3,pD, and (gig) , one finds 

r EW = r s E ^(l + I 2 + 21 cos 6' EW ) 1/2 , (6.8) 
r%w = (^) SM (l + m 2 + 2 m cosQ 1/2 , (6.9) 

T t = : ^2TT- (6-10) 



Since ~ 0.01, the enhancement of r EW remains quite limited in SUSY 

models and it is impossible to enhance it to be of order one. Hence, the contribution 
of r EW to A^£ w+ is negligible respect to the contribution of tew to A^L^. To 
overcome this problem and get the desired relation between A ( ^L n+ and A^^ a 
kind of cancelation between tt and tew contributions to A^L^ is required. Such 
cancelation can be obtained naturally without fine tuning the parameters if ry ~ 
tew-i i- e - the total value of < r^ M . This could happen if k is not very small. 
Therefore, one would expect that the scenarios with dominant chargino contribution, 
where k = 0.061 or k = 0.0019 will not be able to saturate the experimental results 
of A ( ^L it+ and A^ 1 !^ simultaneously. This observation is confirmed in Fig. 3 (top- 
left), where the results of ^4^-^+ is potted versus the results of ^--^o for {k, I, m} = 
{0.061,1.11,2.62} and the other parameters vary as follows: <5j = — 7r,7r/2, and n. 
The angles 6ew an d EW e [—Tt, Tt]- Also dp is assumed to be in the region [tt/4, it/2]. 
Note that in this plot we have taken the Aj£ as constraint. Thus all the points in 
the plot correspond to consistent values of Aj^ with the experimental results. 

Now we consider the second scenario with dominant gluino contribution, i.e., 
(S LR )23 — 0.005e J7r / 3 , {S LL ) 32 = (5^l)32 = 0. In this case, one finds that the maxi- 
mum values of {k, I, m} are give by {k, I, m} = {0.18, 0.014, 0.13}, hence ry is reduced 
from r^ M ~ 0.2 to Tt — 0.12, while tew and r EW approximately remain the same 
as in the SM. In Fig. 3 (top-right) we plot the CP asymmetries AS_ + and Ajf. 
in this scenario, with varying the relevant parameter as before. It is remarkable 
that large number of points of the parameter space can simultaneously accommo- 
date the experimental results of these CP asymmetries. It is slightly surprising 
to get the values of the CP asymmetries A^ within the experimental range, i.e., 
A^L^ e [-0.075,-0.151] and A^L^ G [-0.04,0.12], by just one mass insertion in 
dominant gluino models. This is contrary to the R c — R n results which need gluino 
and chargino combination in order to be within the experimental range. This result 
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Figure 3: CP asymmetry of B — > K ir + versus CP asymmetry of B — > K tt 
for {k, I, m} = {0.061, 1.11, 2.62}, {0.18, 0.014, 0.13}, {0.24, 1.12, 1.48}, {0.32, 0.95, 2.26} re- 
spectively from left to the right and top to bottom. Strong phases 5i = —tt,tt/2, tt and CP 
violating phases 6ew an d &ew res ide between —tt and tt. Finally 6p is assumed to be in 
the region [tt/4, tt/2]. 

can be explained by the cancelation that occurs in A^. w0 between the tt and tew 
contributions and the negligible effect of r^ w to A c K p n+ . 

To be more quantitative, let us consider the following example where (Sf R )23 — 
0.005 e iw ^ 3 and (^ L ) 32 = (^ L ) 32 = 0. In this case one get r T = 0.12, r EW = 0.13 
and r% w = 0.01. Therefore, the main contribution to A^. n+ is due to the linear 
term in tt which is tt sin 5t sin(9 p + 7). With 6p ~ tt/3 and 5t ~ — tt/4, this 
contribution leads to A^Z n+ ~ —0.113. Since r T gives the same contributions to 
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Afj[/L a significant positive contribution from tew is required to change the A^L^o 
and make it positive. With r EW = 0.13, the A^L^ is approximately given by 
A^ p n0 ~ —0.113 + 0.26 siuSew sin(#p — Oew)- It is worth mentioning that although 
9'p and 0' EW are equal in case of single mass insertion, the value of Op and 6ew are 
different due to the different value of k and /. In this example, it turns out that 
6p — Oew ~ tt/9. Hence, one gets A^L^ ~ —0.113 + 0.22 sin Sew- So that for 
Sew ~ 7r /4, one finds A^L ~ 0.04 which is the central value of the experimental 
measurements reported in Table 1. 

We turn to the contributions from two mass insertions: {S LR )23 and {Spp)32, 
which reflect simultaneous contributions from the penguin diagrams with chargino 
and gluino in the loop. Applying the b — > 37 constraints on these mass insertions, 
the maximum values of {k,l,m} is found to be {0.24, 1.12, 1.148}. In this case we 
obtain rp = 0.11, r E w = 0.54 and r EW = 0.06. Therefore, the CP asymmetry Aj^L^ 
is dominated by r E w contribution and in order to get A^L^ of order 0(0.04), a 
small value of the strong phase 5 E w should be used. This makes the possibility of 
saturating the results of A^L^ + and A^-^ is less possible than the previous case. In 
Fig. 3(bottom-left), we present the results of this scenario for the same set of input 
parameters used before. This figure confirms our expectation and it can be easily 
seen that it has less points of the parameter space that account for the experimental 
results of the CP asymmetries than Fig. 3 (top-right). Note also that with two mass 
insertions, the phases 0' P and 0' EW can be considered independent, hence the angles 
Op and Oew are also independent. 

Finally we consider the case of three non- vanishing mass insertions: (5^)23, 
(5^)32 and (S LL )s2- Including the b — > 57 constraints we find that maximum value 
of {k,l,m} is found to be {0.32,0.95,2.26}. The corresponding values of r» are 
r T = 0.10, r E w — 0.48 and r EW = 0.09. It is clear that r T and r E w are slightly 
changed than the previous scenario, while r EW is enhanced a bit. In this case, it will 
be easier to accommodate for A^Jl n0 . The numerical results for this scenario are given 
in Fig. 3(bottom-right) for the same set of parameter space used in previous cases. 
As can be seen from this figure, the probability of accommodating the experimental 
results of different CP asymmetries in this class of models is higher than it in models 
with two mass insertions. However, it remains that the model with dominated gluino 
contributions provides the largest possibility of saturating the experimental results 
of CP asymmetries of B — > Kir. 

6.2 SUSY contributions to the mixing CP asymmetry in B — > i^°7r° 

We turn our attention, now, to the mixing CP asymmetry of B — > K°7i°. As men- 
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tioned before, this decay is dominated by b — > s penguin. Thus, within the SM, the 
CP asymmetry S^o^o should be very close to the value of sin 2(3 ~ 0.73. However, 
the current experimental measurements summarized in Table 1 show that Sk°t° is 
lower than the expected value of sin 2(3, namely 

S K o w o ~ 0.34 ±0.28. (6.11) 

In this section we aim to interpret this discrepancy in terms of supersymmetry con- 
tributions. It is useful to parameterize the SUSY effects by introducing the ratio of 
SM and SUSY amplitudes as follows 



= R n e i6 ^e iS % (6.12) 

Kit 



A svsy (B^K«tt°) 
A SM (B->K°iT ) 



and the angle 9 n is the 



where R w stands for the absolute value of 

SUSY CP violating phase. The strong (CP conserving) phase 5 n is defined by 8^ 
fiSM _ $susy _ rj^- s p arame ^ r i za ^i on j s analogously for those of Sk^> and Sk v ' [5, 
7]. Using this parametrization, one finds that the mixing CP asymmetry Sko„o in 
Eq. fl6.2Q takes the following form 

sin 2(3 + 2R W cos 8 n sm(9 w + 2(3) + R\ sin(2^ + 2(3) 

^>K -k = o ' (6.13) 

1 + 2R W cos 8 W cos 6 W + R w 

Assuming that the SUSY contribution to the amplitude is smaller than the SM 
one i.e. R^ <C 1, one can simplify the above expressions as: 



S K07r o = sin 2/3 + 2 cos 2(3 sin 6 n cos 8 n R w + 0(R 2 n ) . (6.14) 

In order to reduce Sko^ smaller than sin 2(3, the relative sign of sin 6^ and cos 5^ has 
to be negative. If one assumes that sin cos 5„- ~ —1, then R n > 0.2 is required in 
order to get S^o^o within la of the experimental range. 

In the QCDF approach, the decay amplitude of B — > K°ir is given by Eq. (|2.8|) . 
As in the case of B — > <fi(rf)K [5], we will provide the numerical parametrization of 
this amplitude in terms of the Wilson coefficients Ci and Ci defined according to the 
parametrization of the effective Hamiltonian in Eq. fl2.1|) 

* i 

where the operators basis Qi and Qi are the same ones of Eq. fl2.lp . By fixing the 
hadronic parameters with their center values as in Table 1 of Ref. [11], we obtain 

A{B - K\°) = -i^ m \Fl E Hi(ir)(Ci - Q), (6.16) 
^ i=i.. 10,77,89 
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where 



#i(tt 
H 2 (n 
H 3 (tt 
H 4 (ir 
H 5 (n 
# 6 (tt 
H 7 (n 
H 8 (n 
H 9 (ir 



H, 



7T 



-0.7 + 0.00032, 

-0.21 + 0.037i - 0.006X^, 

0.22 - 0.076i + 0.0045X A + 0.0003X A + 0.0065X H , 

0.68 - 0.078i, 

0.2 - 0.001X A + 0.004X A 

0.68 - 0.078i - 0.007X A + 0.014X A , 

0.95 + 0.0004X A - 0.0014X1, 

-0.068 + 0.08i + 0.002X A - 0.0047X A - 0.009X H , 

-1.16 + 0.026i - 0.0015X A - 0.0001X A - 0.003X H , 

-0.67 + 0.08i- 0.0096X H , 

0.0004, 

-0.045. 



(6.17) 



The different sign between Cj and Cj appearing in Eq.( |6.16 ) is due to the fact that 
{K°ir \Qi\B) = —(K°TT°\Qi\B), since the initial and the final states have different 
parity. Comparing the coefficients H^tt) with Hi((p) and Hi(rf) in [5], one finds that 
the Wilson coefficients in these decay amplitudes are different. Thus it is naturally 
to have different CP asymmetries S'^o 7r o, Skq and Sk v ', unlike the SM prediction. 

In order to understand the dominant SUSY contribution to the CP asymmetry 
Sko-kO, it is useful to present a numerical parametrization of the ratio of the amplitude 
R-x in terms of the relevant mass insertions. For the usual SUSY configurations that 
we have used in the previous sections, we obtain 



R n ~ {0.02 x e- MA (6 d LL ) 23 - 40.4 x e~* a01 {6 d LR ) 23 } - {L ^ R} 

- 0.08 x e-* - 013 (5« x ) 32 . 



+ 0.15 x e 



(6.18) 



From this result, it is clear that the largest SUSY effect is provided by the gluino 
contribution to the chromo-magnetic operator which is proportional to (Slr) 2 s anc ^ 
(Srl) 23 - For (5lr) 23 — 0.006 x e 47r//3 and all the other mass insertions set to zero, 
one finds Sko^o ~ 0.34 which coincides with the central value of the experimental 
results reported in Table 1. It is important to note that with such value of (^lr) 23 
the gluino contribution can account for the CP asymmetries Sk4> and Skt, 1 as well [5]. 
Furthermore, if we consider the scenario where both chargino and gluino exchanges 
are contributed simultaneously, the result of R n is enhanced and we can get smaller 
values of Sko^o. 
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7. Conclusions 



In this paper we have analyzed the supersymmetric contributions to the direct and 
mixing CP asymmetries and also to the branching ratios of the B — > Kir decays in 
a model independent way. 

We have shown that, in the SM, the R c — R n puzzle which reflects the discrepancy 
between the experimental measurements of the branching ratios and their expected 
results can not be resolved. Also the direct CP asymmetries A^o^- and A^^o are 
very small while A^o^o and A^L + are of the same order and can be larger. These cor- 
relations among the CP asymmetries are inconsistent with the recent measurements. 
Moreover the mixing CP asymmetry Sk<> w o, which is expected to be sin 2/3, differs 
from the corresponding experimental data. The confirmation of these discrepancies 
will be a clear signal for new physics beyond the SM. 

We have emphasized that the Z-penguin diagram with chargino in the loop and 
the chargino electromagnetic penguin can enhance the contribution of the electroweak 
penguin to B — > Kn which is supposed to play a crucial role in explaining the above 
mentioned discrepancies. We, however, found that these contributions alone are not 
enough to solve the R c — R n puzzle. It turns out that a combination of gluino and 
chargino contributions is necessary to account for the results of R c and R n within the 
b — > 57 constraints. Nevertheless, our numerical results confirmed that the general 
trend of SUSY models favors that the experimental result of R c goes down. 

We have also provided a systematic study of the SUSY contributions to the direct 
CP asymmetries for B — > Kir decays. We found that a large gluino contribution is 
essential to explain the recent experimental data. It is worth mentioning that a large 
gluino contribution is also important to accommodate another controversial results 
measured in the B factories, namely the mixing CP asymmetries S^k and S v 'k- 
Unlike the R c — R n puzzle, we found that the CP asymmetries A^F can be saturated 
by a single mass insertion (5f jR )23 contribution. It has been noticed that a large 
electroweak penguin is less favored by the CP asymmetries A^f- Therefore, one needs 
to optimize the gluino and chargino contributions in order to satisfy simultaneously 
the branching ratios and the CP asymmetries of B — > Kir. 

Finally we have considered the mixing CP asymmetry S K o w o. We found, as in 
S^k and S v 'k, that the gluino contribution through the LR or RL mass insertion 
gives the largest contribution to S^o^o. On the other hand, it is quite possible for 
the gluino exchanges to account for S^o^o, S^k and S v 'k at the same time. 
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